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Abstract. The hydrophobic effect is the dominant force which drives a protein 
towards its native state, but its physics has not been thoroughly understood yet. 
We introduce an exactly solvable model of the solvation of non-polar molecules 
in water, which shows that the reduced number of allowed configurations of water 
molecules when the solute is present is enough to give rise to hydrophobic behaviour. 
We apply our model to a non-polar homopolymer in aqueous solution, obtaining a 
clear evidence of both "cold" and "warm" collapse transitions that recall those of 
proteins. Finally we show how the model can be adapted to describe the solvation 
of aromatic and polar molecules. 
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1. Introduction 

The hydrophobic effect, namely, the free energy cost that non-polar 
solutes pay when transferred into water, is believed to be the dominant 
driving force of protein folding In the native state of real proteins, 
in fact, non-polar residues are buried in the interior of the structure, 
thus minimizing the exposure to water and the subsequent free energy 
cost. 

Nonetheless, the physical properties of liquid water underlying the 
phenomenon of hydrophobicity and giving rise to the characteristic 
behaviour of the change of the thermodynamic functions upon solva- 
tion of a non-polar compound (not only the free energy increase, but 
also the characteristic temperature-dependence of the excess specific 
heat, of the internal energy and of the entropy) are not yet completely 
understood, despite the extensive studies devoted to their investiga- 
tion, ranging from simplified models to numerical simulations (see 
Ref . Q and references quoted therein) . In many studies hydrophobicity 
has been related to the ordering of water molecules around the solute 
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1^, H], but the question is still controversial, since this behavior, indeed 
detected in simulations could be a by-product and not the origin of 
hydrophobicity. An alternative explanation suggests that hydrophobic 
behaviour is related to the process of opening a cavity in water to insert 
the nonpolar solute, and to the interaction between water and solute 0|. 
The decrease in entropy would not be related to bond-induced ordering 
of water molecules, but to the opening of the cavity. 

Here we discuss a model of hydrophobic solvation which is able to 
test the "ordering" hypothesis directly. First, we show that the smaller 
number of allowed configurations of water molecules around a non- 
polar compound is enough to produce hydrophobic behaviour. Then 
we apply our model, which we keep as simple as to be analytically 
integrable, to the case of polymer solvation, and recover for a nonpolar 
homopolymer both "cold" and "warm" collapse transitions that recall 
those of proteins 0, thus strengthening the idea, already put forward 
in Refs. |8|, |9|, that an explicit, though simplified, description of water 
molecules around non-polar solutes can provide a framework for a uni- 
fied treatment of both the "warm" and the "cold" collapse transitions 
of polymers and proteins in aqueous solution. We recall that most 
treatments of protein folding, where water is not taken into account 
explicitly and hydrophobicity is described through effective potentials, 
do not allow a description of cold unfolding. Thus the model introduced 
here might well be relevant to the problem of protein folding. Finally we 
show how our model can be adapted to describe not only the solvation 
of purely non-polar (aliphatic) molecules, but also of aromatic and 
polar ones. This opens the road towards a definition of a model of the 
solvation of a real protein. 



2. Modeling hydration water 

Let us now introduce our model. Consider the water molecules belong- 
ing to the hydration shell around an isolated solute: the solute will 
affect their geometric arrangement, causing a stronger spatial correla- 
tion than in the bulk case. We describe water at the scale of the cluster 
of molecules which are spatially correlated in the presence of a nonpolar 
solute: interactions between clusters are not explicitly considered and 
correlations are lost above this scale. This is consistent with the ob- 
servation ||7| that the contributions to thermodynamic functions from 
the different chemical groups of a residue's side chain are to a good 
extent additive, suggesting the existence of a length scale up to which 
molecules are spatially correlated due to the solute, while correlation 
is small and can be neglected on larger scales. 
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Figure 1. Schematic picture of a hydration water cluster according to the model 
described in the text. The cluster molecules form bonds that can be divided into 
external (e), on-shell (s) and internal (i) bonds. In the external region molecules of 
other clusters can be found, while the internal region contains either a solute (shell 
case) or other water molecules (bulk case). Bonds can be broken (with probability 
pi: see text) if their energy is larger than a given threshold. 

For the sake of simplicity we describe water molecules as two-dimen- 
sional objects with three hydrogen-bonding arms ||l^, |5|, representing 
the projection on a plane of the tetrahedral coordination of a real water 
molecule. The three arms are equivalent: no distinction is done between 
hydrogen donors and acceptors. A schematic picture of a hydration 
water shell according to our model is presented in Fig. |l|. Given a cluster 
of m molecules, we assume that its ground state is characterized by a 
completely formed hydrogen-bond network, not only in the bulk case, 
but also when a solute is present (this involves geometric conditions 
on the solute's shape and size, which we take for granted). Thus each 
molecule in the ground state has 3/2 hydrogen bonds, but their energy 
can be different in the "bulk" and "shell" cases: namely, the energy 
difference per molecule, normalized to a bulk hydrogen-bond energy 
hh, is: 



where hg is the bond energy for shell molecules, and J takes into 
account all the contributions not related to hydrogen bonds (Van der 
Waals and so on). Here h, are positive quantities, while J and K can 
be positive or negative. 
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Our goal is to evaluate the partition function for the cluster in both 
cases: 

Zf-= rd£g.(e)e-^«-'"(^) (2) 

where • = 6, s in the bulk and in the shell case, respectively, is 
the density of states of the cluster at the energy E = hf,e above its 
ground state and Ti^^ is 

n'J'' = hbie + Km6,^s) ■ (3) 

In the above framework, all the important features determining the 
system thermodynamics are encoded in the density of states ^.(e) 
(rather than in the Hamiltonian, which has been substituted by its 
value) . To evaluate g, (e) , we assume that each bond in the system can 
be broken or formed independently of the others, and that every formed 
one can be described by a harmonic potential. A bond is broken if its 
energy[| exceeds h,, while to form a bond one needs also to have a bond- 
ing partner whose arms are correctly oriented. In other words, all the 
configurations with a given number of broken bonds are degenerate in 
energy, and this degeneracy depends on geometrical constraints, which 
can be different in the bulk and shell cases. Thus, at a given energy 
e there will be a probability p,(A,e) to break A bonds, and a given 
number of broken bonds can be obtained with a set of configurations 
of water molecules whose number depends on the presence or absence 
of a nonbonding solute (see Fig. ^). The important point is that when 
few bonds are broken the available number of configurations is smaller 
in the shell case rather in the bulk one, i.e., the solute indeed forces 
an "ordering" of the water cluster. We shall see in the following that 
this fact, together with a different value for hs and hf,, is enough to get 
hydrophobic behaviour. 

Hence, a reasonable ansatz for the functional form of g,{£) appears 
to be: 

/3 \ 

A=0 ^ 

where uj{\,£) is the density of states of the system of harmonic os- 
cillators resulting when A bonds are broken, and 7, (ly) is the number 
of geometric arrangements of the water molecules allowing v unbroken 
bonds. Assuming independence of the bonds, we can write: 

P.(A,e) = fHp^(l-Pi)^"^-' (5) 



^ Considering independent bonds allows one to speak of the "energy of a bond" , 
even if in principle one can only speak of a mean energy per molecule. 
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Figure 2. An example showing the dependence of the "geometric" contribution to 
the density of states <?,, 7, in Eq. on the presence of the solute. To keep the 
correct number of bonds per molecule we fictitiously split water molecules in the 
cluster in two groups; molecules of group a will be able to form only one bond (one 
black arm), while two bonds will be attributed to group b ones. Assuming that bonds 
can be broken independently, in the situation (a) depicted on the left, where both 
molecules are in state 1 (with two on-shell and one external arm; see text), every 
bond may be broken or formed. If we want to evaluate the number of configurations 
allowing v intact bonds, this particular arrangement will be counted for every value 
of v ranging from to |m (three, in this case). But if we rotate one of the molecules 
(right picture) to state 2 (with one internal arm), the shell bond is always broken 
for geometric reasons. The missing bond is definitely lost in the shell case (6), but 
might be recovered by the internal arm in the bulk water case (c) This means that 
this configuration will contribute to the 1/ = |m case just in bulk water. When the 
solute is present, at least one bond is broken, so this configuration will contribute to 
the V = |m — 1 case (1 bond broken) down to the all-broken case = 0, but not to 
the completely bonded v — |m. This example helps one to understand why 7^ < -yt 
when few bonds are broken. 

where pi is the probabihty that a bond acquires an energy larger than 
/i, and breaks, when the cluster energy /i^e is equipartitioned on D 
degrees of freedom, so that pi = e~^'^'^ and T = Notice that 

in principle D = D{m,X) also depends on A, as well as cj(A,e) does: 
the density of states of interacting molecules depends on how many 
molecules are bonded and how many are free, which in turn is a function 
of A. Yet, for the sake of simplicity, we assume the degrees of freedom 
to be always those of bonded molecules, so that D = 2mf and 

uj{X,e) ^iu{e) =Ce^-'^ , (6) 

where / are the degrees of freedom of one molecule (/ = 3 for the 
2-dimensional case), and C is a constant. 
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To estimate 7, (z^) in the 2-dimensional case we consider a clus- 
ter and divide the space in three regions: the cluster's molecules are 
aligned in the shell region, that separates the external one, where water 
molecules of other clusters are found, from the internal one, which 
contains either the nonbonding solute or other water molecules in the 
"s" or "b" cases, respectively. Thus, we talk about "external", "on- 
shell" and "internal" bonds, according to the region they point towards 
(see Fig. Moreover, we classify the orientation of a molecule of the 
cluster according to the direction of its arms, considering, for the sake 
of simplicity, only two states: state 1, with two on-shell arms and the 
third pointing outside, and state 2, with two external and one internal 
arms. In the latter, the internal bond will always be broken when the 
solute is present, while in the bulk water case it can be formed with 
a probability accounting for geometric restrictions on the orientation 
of the internal water. We assume no restriction on external bonds: 
they can always be formed. In this framework, 7, (i^) is related to the 
number of arrangements of the m molecules in the two states allowing 
for 1^ bonds to be formed, and in the end one gets: 



Us=0 Uc=0 



(7) 



where i^i = u — Ug — u^, 



s=0 k=0 3=0 i=0 



(8) 



with /cmax = niin(s, m — s — i^s) and 

In the above equation the following definitions hold: 



7r(s, k) = — 



7rh(s;n, 5, A^) 




(10a) 



(10b) 



(10c) 
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Let us now explain how these expressions can be derived. To keep 
the right number of bonds per molecule (i.e., 3/2), let us fictitiously 
split the m molecules into two groups of m/2 each. Group a-molecules 
will be able to form one bond: a shell one when in state 1 and an 
external one when in state 2. Group 5-molecules will be able to form 
two bonds: a shell one and the external one in state 1, an external and 
the internal one in state 2 (See. Fig. |2|). Hence, state 1 of both groups is 
able to form shell bonds, while internal bonds are possible only in the 
(b, 2) state. Given the number s of molecules in the state 2 (not shell- 
bonding) and the number k of groups of state 1-molecules between state 
2-ones, the probability of making x = itt- ~ s — k shell bonds is tt{s, k) 



given in Eq. ( 10a ). The probability that one of these configurations 
also has i internal bonds depends, first of all, on the probability of 
fishing out j (6, 2)-molecules among the total of s in state 2, given the 
total number of molecules m and the total number of 6-molecules m/2, 
so that one gets a hypergeometric probability - defined in Eq. ( |10bD 
- 7rh(j;s, Y,m). This is not enough, though, because, for geometrical 
reasons, some of the internal arms will not find a bonding partner. 
Indeed, we assume that there are just q, (out of m) positions where 
internal bonds may actually be formed. Their number distinguishes 
the bulk from the shell-water case: Qs = with a non-bonding solute, 
while in the bulk case < qb < m. Again we have a hypergeometric 
probability of placing i of the j molecules with internal arms, in the 
q, good positions for bond formations, on a total of m possibilities: 
7rh(i; j, q,,m). The product of the above probabilities gives the fraction 
of the total number of conformations that is able to form x shell bonds, 
i internal ones and m/2 + s — j external ones, assuming that all the 
molecules with external arms - i.e., group b and state (a, 2) - form 
external bonds. If we now let the bonds be also broken, any geometric 
arrangement of the m molecules allowing v bonds also contributes to 
the cases where fewer bonds are formed. Indeed, given a geometric 
arrangement, we can choose to keep or break the external, internal 
and shell bonds with the binomial probability vTb defined in Eq. ( |lOcD , 
whence the expression of ^, in Eq. (^). Upon summing over all the 
geometric arrangements that can contribute to a pattern with bonds, 
we obtain the expression (0) for 7, {u). 

Upon substituting Eqs. (^J^,0) into Eq. (|^) and then into Eq. 
we get: 

where 

A.=Cy:C'T]7.(lrn-X]l.{\), (12) 
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m-lV 3m/2-A /3rn_\\ 

^- = T^^^.o + E r ■ (-lF2af/^^,,(r.), (13) 
where Kn{z) is the Bessel--fC function and the following definitions hold: 



h.D (j + A) 

(3hl 



^» = ^r2 ' (14) 



r. = 2J(3h,D (i + A) . (15) 



3. Hydrophobic solvation 



Having defined a model for hydration water which should take into 
account the main physical ingredients relevant to the thermodynamics 
of solvation, let us now apply our model to two physically interesting 
cases which can be dealt with in a purely analytical way. To start with, 
we consider the solvation of a single non-polar compound, then we 
apply our model to a non-polar homopolymer. 



3.1. Transfer of a nonpolar solute into water 



Let us consider the cluster of ra water molecules in the bulk case and 
substitute the internal water with a nonpolar solute. This describes the 
transfer of a nonpolar molecule from its gas (non-interacting) phase to 
water, i.e., hydrophobic solvation. We study the difference in thermo- 
dynamic functions. The free energy change is: 

AF = Fsolution - i^water = --^ log X , (16) 

where 

x = ^ e-^^'^'™ . (17) 

Similar expressions hold for energy, entropy and specific heat changes. 
The temperature dependence of these functions, reported in Fig. 
shows the hallmarks of hydrophobic behaviour, even for K = (i.e., 
without any ground-state energy difference). In fact, we find a maxi- 
mum in the free energy cost, a pronounced and positive peak in the 
specific heat difference AC, and minima in both AE and TAS; then, as 
T grows, AE and TAS" cross the zero - thus defining the characteristic 
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Figure 3. Free energy, energy, entropy and specific heat changes upon solvation of a 
nonpolar solute: AF (dash-dotted line), AE, (dashed), TAS (dotted), AC/3 (solid 
line; 1/3 is for rendering purposes). Here m — 4, qt = m/2 (half of the molecules 
can form internal bonds, in the bulk case), hs/hb = 1.2, J = 0.3 (so that the 
ground-state energy shift K = 0). Energies and temperatures are expressed in bulk 
hydrogen bonds units (hb); Boltzmann constant fcs is set to 1 and specific heat 
is, accordingly, adimensional. The shaded region can be compared with numerical 
simulations (e.g., with Fig. 5 of Ref. 0) and experimental results like those reported 



temperatures commonly referred to as Th and Ts - and eventually be- 
come positive. Notice that the peak in AC cannot be recovered within 
the Information Theory Approximation, usually applied together with 



3.2. Polymer in solution 

Let us now turn to the study of a nonpolar homopolymer in solution, 
taking into account just the behaviour of water clusters in the vicinity 
of a monomer, and disregarding interactions between monomers and 
between water clusters. Our goal is, in fact, to understand the effect of 
the hydrophobic effect alone on polymer behaviour. 

We model a polymer as a A^-step self-avoiding walk on a two di- 
mensional lattice with coordination z. On each lattice site there can 
be either a monomer or z clusters of m water molecules, so that each 
monomer-water contact involves one cluster. The Hamiltonian follows 



where N\y = (z — 2)A^ -|- 2 is the highest number of water sites that 
can be in contact with the polymer, and Ij is the number of contacts 



in Ref. §. 



the cavity approach pT[ |. 



from Eq. (^: 




(18) 
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between the j-th water site and the monomers. The partition function 
of the polymer in solution reads as: 

Z = Y,Z{C) = Y,C{n,)Z{n,) (19) 

C ric 

where C are the conformations of the polymer and Z{C) the restricted 
partition function, obtained tracing over water variables at fixed con- 



formation C. Due to the form of the Hamiltonian given in Eq. (18), 
Z[C) depends only on the total number of water-monomer contacts 
ric, and C, {ric) is the number of SAWs characterized by the same value 
of ric- Z can be factorized as 

Z = Z^Zi , (20) 

where Zi, = (^z^^Y^^ = A^^"^ is the contribution of all water sites 
when in contact to other water, and 

2/ = ECK)x"^. (21) 

In the following we shall study the specific heat of the system, that, ac- 
cording to the above factorization, is the sum of a bulk contribution Cb 
and of an interaction contribution C/. We shall also study the average 
number of water- monomer contacts: 

{nc) = x^\ogZi, (22) 

which is a measure of the compactness of the polymer. To exactly 
evaluate the above quantities, an exhaustive enumeration of the SAWs 
should be performed, in order to obtain C, {ric)- However, if we restrict 



ourselves to a square lattice, the numerical estimates reported in |12| 
allow us to write 

C (uc) ^ Co^^(«oiV)"'("^^ exp(aoA^) , (23) 

where w{nc) = — nc)/2 is the number of monomer-monomer con- 
tacts, Co is the number of SAWs of length and ao = 0.164. Such 



an estimate is expected to be very good if A^ is large |12|. Hence, an 
analytical expression can be found for Zf. 

Z, = ^x^-r(u;_ + 1, ^)e--^(i"^) , (24) 

where r(n, x) is the incomplete F-function and t^max = (A^VF — n™™)/2 
(n™™ accounts for the fact that the smallest number of water-monomer 
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Figure 4- Excess specific heat Ci/Nw of tfie nonpolar fiomopolymer in solution, for 
different lengtfis: JV = 2 x 10^ (dastied), iV = 10** (dasti-dotted), iV = 10* (solid 
line). Parameters as before. 




0.2 0.4 0.6 



Figure 5. Normalized average number of water-monomer contacts, fc = (nc)/Nw, 
for different polymer lengths (the same as before). Parameters and rendering 
conventions as before. 



contacts is determined by the length of the perimeter of the globule: 

The results obtained C/ and (ric) stemming from Eq. (^) are re- 
ported in Figs. 1^ and ^, for different lengths of the polymer. 

The presence of both "cold" and "warm" collapse transitions, sig- 
nalled by the drop of the number of contacts and by the jumps in the 
specific heat is strikingly evident from our results. As N grows both 
transitions get sharper, thus suggesting the existence of a true phase 
transition in the thermodynamic limit. This phenomenology is very 
close to that of proteins in solution: our results confirm that, as already 
suggested in Refs. ^, the explicit treatment of the solvent - though 
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in an extremely simplified way - is a natural and powerful way to get a 
unified modeling of both the cold and the warm unfolding transitions in 
proteins. Moreover, it is worth noticing that a cold collapse qualitatively 
very similar to that predicted by our model has been experimentally 
observed also in homopolymers like the poly(N-isopropylacrylamide) 



4. Solvation of aromatic and polar compounds 

Proteins are heteropolymers whose monomers are aminoacids. The 
aminoacid residues (i.e., those parts of the aminoacids which can be ex- 
posed to water once they are linked with other aminoacids in a protein 
chain) can be roughly split into two classes: polar (or even charged) 
and non-polar. However, as far as the qualitative features of their 
interaction with water is concerned, not all the non-polar molecules 
behave in the same way. One can identify two families, aliphatic and 
aromatic residues |^, with different thermodynamic behavior. Aliphatic 
compounds are "purely" hydrophobic, i.e., show a positive jump in both 
the specific heat and the free energy upon solvation. On the contrary, 
when an aromatic residue is transferred into water from the gas phase 
there is a positive jump in the specific heat but there is a free energy 
gain too, i.e., AF < 0. Aromatic molecules show, in a way, an interme- 
diate behaviour betweeen hydrophobic and hydrophilic. Moreover, also 
in the case of polar molecules there is a change in the thermodynamic 
functions: AF < and AC < 0, so that they can be reasonably called 
hydrophilic. 

The model of solvation discussed so far works very well when ap- 
plied to purely hydrophobic solvation, i.e., to the solvation of aliphatic 
monomers or homopolymers. However, in view of a generalization of the 
approach presented here to heteropolymers and especially to proteins, 
one should devise a way to model also the solvation of aromatic and 
polar residues. As we are going to show in the following, our model 
is able to accomplish both tasks, provided one chooses the parameters 
according to the physical requirements of these situations. This is not 
only very promising for future applications to proteins, but also tells 
us that our model grasps, though in a very simplified way, some of the 
basic physics of water-solute interaction. 

4.1. Aromatic molecules 

Figure ^ presents trends in the changes in thermodynamic functions 
upon solvation that qualitatively recall those obtained experimentally 
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Figure 6. Free energy, energy, entropy and specific heat changes upon solvation of 
an aromatic solute: rendering conventions as in Fig. ^. Here m = A, qb = rn/2, 
hs/hb = 0.98, J = —0.07. The shaded region can be compared with experiments 
reported in Ref. Q. 



for aromatic residues |^. The plots are obtained with hs/hb = 0.98, 
J = —0.07, Qs = 0, Qb = 772-/2; more in general, hs/hb ^ 1 and Qb = m/4 
would provide similar results, as far as the experimental temperature 
window is concerned (data not shown). In the framework of the model, 
this implies that hydrogen bonds for bulk and hydration water are 
more or less equivalent, and that shell molecules are less likely to form 
internal bonds than in the bulk case (or cannot at all). At the same time 
there is a small gain, related to J, in hydrating the solute, suggesting 
that other interactions than hydrogen bond could favour the aromatic 
hydration. A clear intepretation of this fact is not straightforward: it 
is unlikely that a negative J could come from Van der Waals contri- 
butions prevailing over the electrostatics of bulk water: electrostatics 
interactions between solute and solvent should come into play, instead. 
It is difficult to say if their effect could be described by an isotropic 
term J; we cannot exclude also that a negative J compensates the fact 
that in the model internal bonds have the same energy as external 
and on-shell. Even if the precise scenario could be determined only 
by all-atoms simulation, the hypothesis that electrostatic interactions 
between solute and water are more important for aromatic than for 
aliphatic compounds is supported by the experimental evidence that 
aromatic residues present a partial polar nature, with an excess negative 
charge on the faces of the aromatic ring and a partial positive charge 
of the hydrogen atoms at the border ||^, |l^, ^. 
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4.2. Polar molecules 



Figure |^ presents thermodynamic trends which are in a good qual- 
itative agreement with experimental findings for polar residues 0. 
Again, the choice of the parameters is not unique, since the region 
to be compared to experiments is not as sensitive as the low T one to 
parameters changes, and a comparison of a two-dimensional model to 
real experiments cannot be made truly quantitative. Nevertheless, there 
are some crucial aspects characterizing polar behavior: a low value of 
hs/hf, (i.e. less than 0.6) is required to produce the correct trend in 
AC together with a low value for T^, whence the relation AE < AF. 
On the other hand, a low ratio hs/hf, would yield a positive energy 
shift K] thus, it must be compensated by a large and negative J, 
favouring solvation, if the experimental values for a large and negative 
AF are to be recovered. It is important to notice that, for low values 
of hs/hf), the different degeneracies 56 (e), 5s (e) in bulk and shell water 
come essentially from the different probabilities of breaking bonds: the 
geometric arrangement of the molecules almost plays no role. 

The need for a low ratio hs/hb and J < suggests a possible phys- 
ical explanation of the parameters: the presence of a polar compound 
hinders the formation of hydrogen bonds between water molecules by 
competing with it; a strong direct solute-water interaction (J is almost 
equal to hf,) results in weakened hydrogen bonds among shell water 
molecules. 
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5. Concluding remarks and future developments 

Wc have introduced and diseussed a statistical-mechanieal model of 
non-polar solvation in water. In our model the water degrees of free- 
dom are explicitly considered, though in a very simplified way: water 
molecules are modeled as planar objects with three equivalent bonding 
arms (roughly speaking, the projection of water's tetrahedral coordina- 
tion into a plane) and all the treatment is for a two-dimensional system. 
Nonetheless, this allows us to obtain an exact analytic solution of the 
model, which unambigously shows that hydrophobic behaviour can be 
related to a combination of the strengthening of hydrogen bonds to- 
gether with the reduction of available configurations for water molecules 
when the solute is present. Our model thus supports the conjecture 
that the physical origin of hydrophobicity is the formation of "ordered" 
cages of water around the solute, with stronger hydrogen bonds than 
in bulk, even if we cannot say which one of the two ingredients above is 
more determinant. Applying our model to a hydrophobic homopolymer 
we have clearly shown that it is swollen at low temperatures, then 
becomes maximally compact and eventually unfolds again at higher 
temperatures. This behaviour closely reminds that of proteins, which 
undergo both a cold and a warm unfolding transition. 

In view of the application of our model to proteins, we have shown 
that it can be adapted to describe also the solvation of aromatic and 
polar molecules. To study the behaviour of a model of a real protein 
along these lines, however, it will be necessary to extend the model to 
the three-dimensional case and to determine, by means of a comparison 
with experiments, the right parameters for all the aminoacids. Work is 
in progress in these directions. 
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